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Abstract 

In three dimensions there is a logarithmically divergent contribution to the entanglement entropy 
which is due to the vertices located at the boundary of the region considered. In this work we find 
the corresponding universal coefficient for a free Dirac field, and extend a previous work in which the 
scalar case was treated. The problem is equivalent to find the conformal anomaly in three dimensional 
space where multiplicative boundary conditions for the field are imposed on a plane angular sector. As 
an intermediate step of the calculation we compute the trace of the Green function of a massive Dirac 
field in a two dimensional sphere with boundary conditions imposed on a segment of a great circle. 

1 Introduction 

A typical manifestation of entanglement is the presence of correlations for the results of measurements 
on commuting sets of observables. In quantum field theory, the vacuum fluctuations provide a source of 
quantum entanglement between spatially separated regions. A natural measure of this entanglement is the 
geometric entropy, which is defined as the von Neumann entropy S{V) = — tr(/9y log /3y) of the reduced 
density matrix of the vacuum state pv = tr_y |0) (0| corresponding to V. The trace in the expression for 
pv is taken over the Hilbert space generated by the degrees of freedom lying outside V. 

Recently, there has been an increasing interest in understanding the manifestations of entanglement 
in quantum field theory (QFT). This is in part due to the conjecture that black hole entropy may be- 
come geometric entropy once the role of gravity is elucidated [1]. The investigation has also revealed 
that the entanglement entropy is a very interesting tool in the study of certain aspects of QFT. Very 
different phenomena such as quantum phase transitions |2j, confinement [3j, topological phases [4J, and 
the renormalization group irreversibility [5], leave an imprint on the behavior of the geometric entropy. 

Our hope is that this quantity, which has very powerful non perturbative properties [6], may help 
to classify the different QFT, or even provide an alternative setting where QFT may be expressed in 
a geometrical way. In this sense, it is remarkable that there is a concrete proposal in the context of 
the AdS-CFT duality where the entanglement entropy of the boundary CFT may be computed by purely 
geometric means [7]. This hypothesis has been supported by further arguments ^ and has passed already 
several nontrivial tests [9]. 

With this general motivation, in this paper we study the entanglement entropy associated to a spatial 
polygonal set in (2 + 1) dimensions for a free Dirac field. The structure of divergences the entropy develops 
in the continuum limit is given in three dimensions by the expansion 

S{V) = gi[dV] e-' + go[dV] log(eA) + S^{V) , (1) 

where Sq{V) is a finite term, e is a short distance cutoff, and A is a parameter with dimensions of mass. 
The coefficient of the logarithmic term 50 [t^^] is a universal quantity depending on the set boundary 
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dV [10\ lllj. Because of its ultraviolet origin, for general polygonal sets we expect it to consist in a sum 
of contributions from the individual vertices 



9o[dV] = Y,si^i), (2) 



where the Xi denote the vertex angles. The function s{x) satisfy s{x) = s{2tt — x), which is a consequence 
of the symmetry in the entropies S{V) = S{—V) due to the purity of the vacuum state. 

Here we find the function s{x) = s^ix) for a Dirac field. In a previous work we have found some 
quantities related to the function ss{x) corresponding to a free scalar field [11]. The logarithmically 
divergent term in the entropy due to the vertices has also been found in the holographic proposal by Ryu 
and Takayanagi [12], presumably corresponding to an interacting CFT. We will also compare our results 
with this last work. 

The method usually employed to find these entropies is based on the following representation 

S(F)=limJ5i<i^41. (3) 

n^l (1 — n) 

in terms of the trace of powers of the density matrix. Then py is expressed as 

where Zn is a functional integral on an n-sheeted d dimensional Euclidean space with a conical singularity 
of angle 27rn located at the boundary of the set V (see for example |13j). Here d is the space-time 
dimension, d = 3 in the present work. The replicated space is obtained considering n copies of the d-space 
cut along V, and sewing together the upper side of the cut in the k^^ copy with the lower one of the 
{k + 1)**^ copy, for k = 1, n, and where the copy n + 1 coincides with the first one. 

In general, obtaining these integrals explicitly is a very difficult problem since we have to deal with a 
non trivial manifold resulting from the replication method. Moreover, the entanglement entropy follows 
from Z[n] for integer values of n through an analytic continuation to n = 1. In the case of free fields, a 
great simplification follows by mapping the n— sheeted problem to an equivalent one in which we deal with 
n decoupled free fields 0^ for = 1, n [HKTB]. These fields are multivalued and live on the Euclidean 
d dimensional space with boundary conditions imposed on the d — 1 dimensional set V given by 

= e'^'^l"^(r), feV. (5) 

Here (f)^f^^ and 0^. ^ are the limits of the field as the variable approaches V from each of its two opposite 
sides in d dimensions. In this formulation we have for fermions 

(n-l)/2 

log(trp^)= Yl log^[e'2-^]> (6) 

fc=-(n-l)/2 

and for bosons 

n-l 

log(trp^»=^logZp2.^J ^ 
fc=0 

where Z[X] is the partition function corresponding to a field on a single copy of the Euclidean space, which 
is multiplied by a factor A when the variable crosses V. Although Z [e*^'^'^/"] has been found explicitly 
for several cases of interest |1H [T^ [15] , the evaluation of the entropy in this scenario is still limited by 
the difficulties in getting the analytic continuation of the sum for non-integer n in order to take the limit 
([3]). Only very recently [T6|, this problem has been solved by exploiting the relation between ([6]) and ([7]) 
with the expressions for pv in terms of correlators for free boson and fermion systems [17^ [T8] . The result 
is 

poo 
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for free fermions, and 

for free bosons. Note that the functional integral involved has a boundary condition on the cut V given 
by a real factor (instead of a phase factor as in ^ and ([7])), being a positive factor for fermions and a 
negative one for bosons. The calculation of the geometric entropy for the set V is then effectively reduced 
to the one of the corresponding functional integral Z[X]. In particular, the term logarithmic in the cutoff 
appears in this partition function, with a coefficient proportional to the conformal anomaly induced by 
the boundary conditions ([5]) [19j. 

In order to obtain s^ix) we calculate the Dirac partition function for a set V formed by a plane 
angular sector. We first reduce the problem to a two dimensional one on a sphere with a cut along a 
segment of a great circle, and then relate the Green function on this sphere to the corresponding scalar 
problem studied in [11] . The reported results in this last work concern only the exact values for the trace 
of powers of the local density matrix for a scalar, and only the small angle limit for ssix) was obtained. 
Here we use the formula ([8]) which allows us to give the analytic result for the entropy in terms of a system 
of coupled nonlinear ordinary differential equations. 

The outline of the paper is as follows. In the next Section we show how the three dimensional problem 
in Euclidean space can be reduced to a two dimensional one on a sphere. In the third Section, we show 
the correspondence between the Dirac and scalar Green functions on the cut sphere. Then, in Section 
IV, we use these results and the ones of ref. [IT] to find the logarithmic contribution to the entropy. We 
conclude comparing the coefficients s{x) for the scalar and fermionic cases with the data given by Hirata 
and Takayanagi in [12], corresponding to a conformal theory holographic to gravity in four dimensions. 



2 Vertex contribution to the entropy in three dimensions 

In order to extract the contribution of the vertices to the entropy, the simplest set V to consider is an 
infinite two dimensional plane angular sector of angle x. We take its vertex on the origin of coordinates. 
The associated partition function for a Dirac field ^ in three dimensions is 

Z[e'2-'^] = J P^tp^e-/dr-3>I>t253viy ^ ^^^^ 

where P3 is the Dirac operator in three dimensions given by 

Vs = iYd, + f,), (11) 

and 7* = cr* are the Pauli matrices. The boundary condition for the spinors is 

$+(f) = e*2™^-(f) , rGV. (12) 

Here and are the limit values of the field on each of the sides the two dimensional angular 
sector V has in three dimensions. Since by reflexion symmetry on the plane containing V we have 
2-jg«27raj _ ^[g-*27raj choosc without loss of generality 

a G (0,1/2). (13) 

Note that here we are using a real o, which allow for a direct computation of trp". In order to use the 
formula ([8]) to obtain the entropy we will turn to imaginary a at the end of the calculation. 

Due to the boundary condition, the field is singular as we approach the edges of the plane angular 
sector V or to the vertex point. However, in order to have finite action, we have to impose that it diverges 
at most as 



3 



as r approaches dV (the edges of V), where d is the distance between r and dV. Also we have 

^{r)r^\r\''\ K2>-1, (15) 

as we approach the origin of coordinates. The same conclusion can be reached considering that the 
induced current circulating around dV has to be finite. 

The functional Z is calculated exploiting the relation between the free energy and the Green function 
Gg)(f,r-), 

-- fdr' (vl/t(f)^(f)\ = TrGg) , (16) 



where fi is the field mass. The Euclidean Green function Gd satisfies the equation 

V,G^^\f,f') = 6^{f-r'). (17) 

2.1 Dimensional reduction 

The Dirac operator and the boundary conditions allow the separation of the angular and radial equations 
in polar coordinates. We will use this fact in order to reduce the problem to one in two dimensions. In 
this coordinates P3 writes 

Vs = -+rdr+fi, (18) 
r 

with 

D = (j'de + i^d^) , (19) 

and where the new gamma matrices 7 are 

f = r|^a\ ^<t> = r^a\, f = = -isinOj'^^ . (20) 

OXi OXi OXi 

The operator D satisfies 

Df +fD = 2, D^- YD = An , (21) 

where 

Id d 1 52 

sm B 00 00 sm 9 0(p^ 

is the Laplacian on the sphere. 

Since (-^ + Y^r) is anti-hermitian we can write the eigenvalues of P3 as ;U + i/3, with (3 real. The 
corresponding eigenfunctions, 

i^/3,u = (/i + iP)i^f3,u , (23) 

have the general form 

i^f3,u = f{r)'^u{e, 0) + g{r)f^u{0, 4>) . (24) 
Here ^y{6,(j)) are the normalized eigenfunctions of the angular operator 7*^-0, such that 

fD^, = i^^u, j dVt $t = 5^ ^, . (25) 

Since 7'"L) = (7'"!))^^ the eigenvalues v are real. From eqs. (I2ip it also follows that for each eigenvector $ 
with eigenvalue v of this operator, there is another eigenvector 7''<I> with eigenvalue 2 — v. The spectrum 
of Y D is then symmetric around the point 1, and without loss of generality we can take > 1 in equation 
(121. 

The equation ([^H]) reduces to two coupled ordinary differential equations for the radial functions /(r) 
and g{r) 

fir) 



r 



u + f'{r) = i(3g{r), (26) 



^(2-r.)+5'(r) = i(3f{r). (27) 
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The solutions are given in terms of Bessel functions. Keeping only the solutions for /(r) and g(r) which 
satisfy the condition (|15p we arrive at the normalized eigenfunctions 



i^(3,u = i—^^ — ^ — <P) + V ^ Y'^u{0, (f>) , (28) 



dr j dn V^, ^P(,^, = 5,, - . (29) 
Here we have made use of the orthogonality relation for the Bessel functions 

/•oo 

y/W drrJsirP)Js{rp') = 5i(3-p'), (30) 
Jo 

which is valid for s > —1/2. 

Then, the Green function writes 

/■oo 

G^^\r',e',ip',r,e,^) = Y^ / dp^p,,ir',e',^')ip + i(3)-'i;l^{r,e,^). (31) 

The trace is 

POD 

u>l 



(32) 



These integrals do not converge. However, here we are interested in the dependence of trG^^ on the 
angle x and a, which is finite and calculable. In other words, (I32p just requires a subtraction of a global 
constant. In fact the integrals are divergent because the function x I^{x) K^{x) goes to 1/2 for x — > oo, 
and for any v. Subtracting this global constant inside the integrals leads to 



trGg) = -^-15^(^.-1). (33) 

iy>l 



This sum over the angular modes is also non convergent. We subtract a global constant by conveniently 

(3) 

defining trG^ = on x = vr. At this angle there is no vertex and the logarithmic contribution vanishes. 

In the appendix A we show that there are no eigenvalues v in the interval (1/2,3/2). Since the 
spectrum of values of 1^ is symmetric around 1, we use this fact to rewrite the trace as the following sum 
over all the angular eigenvalues 

'■•Gg' = -^E I" - 11 = -^E I" - V2| = -^trirO - 1/21 . (34) 

We make this shift since the operator on the right hand side is directly related to the Laplacian, 

{fD-l/2f = -A + l/4. (35) 

This will allow us to connect the present problem to a similar one of a scalar field in a sphere previously 
studied in [TT]. 

Note that by dimensional arguments the logarithmic divergent contribution to the partition function 
for an infinite plane angular sector must be proportional to log(e;u). This means by (I16p and ()34p that 
the sought logarithmic contribution to log Z is 

log Zliog = -^tr \fD - 1/2| log(e/x) . (36) 
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The trace in (I36p can be calculated using the integral representation in terms of the resolvent [20] 



1 

ti\fD-l/2\ = / dmmtY{i{fD -1/2) +m)~^ . (37) 

J-oo 

The operator 

V2 = i{fD - 1/2) + m (38) 

is a two dimensional Dirac operator on the sphere where the parameter m plays the role of a mass. In 
the next Section we find the trace of the Green function G^, = , which is required by (I37p . 

3 Green function on a sphere with a cut 

In the dimensionally reduced problem we have to find the trace of the Green function of a Dirac field on 
a two-dimensional sphere. This satisfies 

V^G^^^ = ^ 5\z - z') . (39) 

The following boundary conditions for the spinors are imposed on a segment of a great circle (we choose 
this later on the equator) 

lim ^'(^/2 + e,ip)= e*^^'' lim ^!{t:/2 - e, v?) , S [(^1, ^2] • (40) 

e— >0+ e— »0+ 

We leave the end points of the cut Li = (7r/2,(^i) and L2 = (7r/2,(^2) free, and at the end we will put 
(p2 — = X. The same boundary conditions are satisfied by the field derivatives. 
The Dirac operator (I38p on the sphere can be writter0 

V2 = -f^de + ^^d^-'- + m. (41) 

The gamma matrices on the sphere 7°^ = ^7^7" satisfy 

{7^7^} = 25"^ (42) 

and g^^ = 1, g'^'^ = csc^(0), g^^ = g"^^ = are the components of the metric tensor. Explicitly they are 

ie-^M -1 e-^^coi{9)\ 

-ie'^ y ' V e^'^cot(0) 1 J ■ ^ ^ 

vl = -i^'de + ^%-'-) + m, (44) 

V2Vl = -An + ^ + m\ (45) 



The adjoint operator reads 



and we have from (|35 



^The Dirac operator in curved space is 7^(9;j + F^) + m. The spin connection r,j is such that the covariant derivatives of 
the 7" are zero, D^'y'^ = 9^7" + F'j^^'y" + = 0. Here we can further determine it by imposing that 7'^r^ is diagonal. 

This gives the operator j^dg +7*^9^ — i + m, which differs from (|41|l by a constant term i/2. As noted in the previous Section 
we have chosen to make this shift because in this way our fermionic operator is related to a scalar with a real mass. The 
ordinary Dirac operator is related to a complex mass scalar, and in this case we can not use the results of 111' in a direct form 
(this case could also be treated with some obvious modifications of that work). See [21J for further elaborations regarding 
the dimensional reduction of the Dirac equation in three dimensions in polar coordinates. Note also that the present problem 
has been solved for the hyperbolic plane where the partition function of a Dirac field with the boundary conditions analogous 
to (|40p is expressed in terms of Painleve functions [22] . 
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3.1 Relation between the scalar and Dirac Green functions 

Consider a complex scalar field with mass M on a sphere with a cut where the boundary conditions 

.(2; 



analogous to (HO]) are imposed. The Green function Gg {z^ z') satisfies the following equations 



{- An + M^), of {z,z') = ^5\z-z'), (46) 
hm (z, (^/2 + e, ct>)) = e^^™ ^- ^(2) ^^/^ -e,<P)) , [^1,^2] ■ (47) 

In order to explicitly relate this problem to the fermionic case discussed in the previously, we define an 
auxiliary quantity G as 

G{z,z')=VlGf{z,z'). (48) 

Eq. (I45p implies that 

V,G{z,z') = ^ 5Hz-z'), (49) 
where the scalar and the fermion masses are related by 

= i + (50) 

Thus, the difference 

G^^\z,z')-G{z,z') = Q{z,z') (51) 

satisfies the Dirac, 

VQ{z,z')=0, (52) 

and Helmholtz equation, (— + M^)Q(z, z') = 0, without sources. Therefore it would be identically zero 
if it where bounded. Q{z,z') is however unbounded at the extreme points Li and L2 of the cut. Our 
strategy is to find a function of z which also satisfy the Helmholtz equation without sources and have the 
same singularities as Q{z,z') for z going to Li and L2. Then Q{z,z') must be identical to this function 
by the uniqueness of the solution of the Helmholtz equation for bounded functions. 

The singular behavior of Q{z,z') (as a function of z) is due to two different sources. The first one is 

(2) 

because the scalar function G\, {z, z') is bounded at the cut, but the leading terms go as [(9? — 9?i) +i(0 — 
7r/2)]i-'^ and [((/?- (^i)-z(e-7r/2)]" for z Li, and [((^2 - 92) +i(0 - 7r/2)]i-'^ and [((^93 -99) -i(0-7r/2)]" 
for z ^ L2. Thus, their derivatives in G{z, z') are unbounded at these points. Using the results of previous 
works [mdS]) these singular terms in Q{z, z') are computed straightforwardly in terms of functions Si{z), 
S2{z), and the reflected ones, Si{Rz), S2iRz), where R is the reflexion operator which maps Li in L2 

R{e,<P) = {9,^i + ip2-^). (53) 

Here Si{z) and S2{z) are uniquely defined by two requirements. They satisfy the homogeneous equation 

{-An + M^)Si{z) = 0, {-An + M^)S2{z)=0, (54) 

and diverge as z — > Li as 

Siiz) -^ii^- ^1) - - ^m'^ , S2iz) ~ - ^1) + ^{0 - ^/2)r-' . (55) 

Then, this first source of singularities can be taken into account by the term 

Qi{z,z')= - iA7Ta{l-a)(^^'\L,S2{zrSi{z') + j'\L,S2{RzrSi{Rz')^ 

+ 4^a(l - a) (^^^\L,S2{zrSi{z') - j%,S2{Rzr Si{Rz')^ . (56) 

The second source of singular terms in Q{z,z') is that Gj^ {z,z') has itself singularities according to 
(fTil) . These must be proportional to [(99 — Lpi) + i{9 — 7r/2)]^" on Li and [((/?2 — </')+ i(^ - '^1'^)]^"' on 
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in order to respect the boundary conditions and (I14p (recall that a G (0, 1/2)). Thus, we can always 
write these singular terms in Q{z, z') (as a function of z) as a combination of 5*1(2:) and Si{Rz), 

Q2{z,z') = Si{z)* {zq{z')1 + zi{z')ai + Z2{z')a2 + Z2,{z')a-i) 

+ Si{Rzr (yo(^')l + yi{z)<yi + y2{z)a2 + y^{z')a-i) , (57) 

with unknown coefficient functions Zi{z') and yi{z'). 
Therefore we have 

Q{z,z') = Qi{z,z') + Q2{z,z'). (58) 

Then, the functions Zi{z') and yi{z') can be obtained in terms of the functions Si{z') inserting the above 
expression for Q{z, z') in the Dirac equation (fS^ . and using the differential equations for Si{z) and S2{z) 
given in [11] (eqs. (81-84) of that work). Here we need only the part of the trace of G^, = G + Q which is 
odd in the mass m, since the terms even in m do not contribute to the integral (|37|) . After some algebra 
we find 



trG^n\x,m,a) = 2mtrG^"^ - ^7'"' " ' ' . (59) 



.(2) 



167ra(l - a)m (4/?i Xx cos(a;/2) -hBx sin2(x)) 



odd 



M (4/3i2-52sin2(x)) 



Here tr Gg , b, Pi, Xi and Bi are functions of x which were studied in [TT] . They are given in terms of a 
system of algebraic and ordinary differential which we include in the Appendix B. 



4 Results and conclusions 

Now we can write the formula for the coefficient of the logarithmic contribution to the entropy due to the 
vertices for a Dirac field in three dimensions. The result follows combining ([8]), ()36p and (I37p . We have 

roo -j^ roc 

sd{x) = / dt 27 — T / dmm tvG^ {x,m,—it) . (60) 

Jo 2 sinh (vrt) J_oo 

The relevant part of trGoix, m, —it) is given by (1590 and the formulae at the Appendix B, where we have 
to make the replacement a = —it. The imaginary part cancel in (j59p as it should. 

The analogous result for a complex scalar field follows from [llj with the help of the analytic contin- 
uation described in [16j. We have 

ss{x)= dt 2 / dM M (M^ ti Gf{x,M, -it + 1/2) , (61) 

Jo cosh (vrt) J 1/2 

where we have to replace a = —it + 1/2 in the equations (I72II89P of Appendix B. 

An economic way to numerically integrate the equations (j60p or ()6ip is to expand the functions involved 
in (I73II89P in Taylor series around x = n, and obtain analytically the coefficients using the differential 
equations. Then the above integrals over t and the mass can be done for each coefficient separately. 
With this method we have produced the curves of figure 1, which show sd{x) up to order (x — vr)^^. 
In the picture are also plotted the values of for x = 7r/4, Tr/2 and 3/4 tt obtained by numerical 
simulations in the lattice. They show a perfect accord (around one percent error) with the analytical 
results. These particular values of the angle are the ones for which the coefficient can be calculated 
in absolute terms with very small error on a square lattice of limited size (in the present case it was 
100 X 200 points). The numerical methods consist of evaluating the entropy for a massless Dirac field (see 
[14:\ I17j) for a given shape (square, triangle, etc.) and different overall size A, and then fitting the result 
as 5 = Co + Ci A -|- C_i A^^ -I- C-2 A~^ — sd log(A). It is also possible to evaluate very accurately sd for 
specific combinations of angles using rectangular triangles. We also obtain in this case a perfect accord 
with the analytical results. 

In [7J Ryu and Takayanagi propose a purely geometric method to compute the entanglement entropy 
for certain conformal field theory in the AdS-CFT context. In this way Hirata and Takayanagi find 
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Figure 1: The coefficient sd{x) of the logarithmic term in the entropy for a Dirac fermion computed by 
its Taylor series expansion around x = vr up to order 14 (solid curve). The dashed lines are the Taylor 
expansions of lower order, from order 2 to 12. Also shown are three points of the curve at angles 7r/4, 
7r/2, and 3/47r, obtained by direct numerical evaluation of s^ix) in a lattice (see the text for details). 



a logarithmic contribution to the entropy due to the vertices in three dimensions which presumably 
corresponds to an interacting CFT [12]. The coefficient of the logarithmic term snix) in this holographic 
theory is given in parametric form in terms of a variable go £ (0, cxd) as 



x{go) 



SH{go) 



dz 



{z^ + glW{z^ + gl + l){z^ + 2gl + l) 



N 



dz 1 



+ 9l + l 
+ 25^ + 1 



(62) 
(63) 



In order to compare the functions sd{x)^ ss{x) and sh{x) we have to choose a normalization N for 
sh{x) in (j63]) . Remarkably, the quadratic coefficients of sd{x) and ssix) turn out to be equal. We have 
then chosen to normalize sh{x) such that its quadratic coefficient coincides with the ones of sd{x) and 
ss{x). In Table 1 we show the first four non zero Taylor coefficients of so{x)-, ss{x) and sh{x) for x 
around vr. 

As explained in |llj, all of these functions diverge as 1/x for x going to zero. In the Table 1 we also 
show the coefficient of 1/x in the expansion around x = 0. For sd{x) and ss{x) this small angle behavior 
is given by jTl] 

!^dtC{t)/7: 



s{x) 



X 



(64) 







C4 








s{7t/2) 


s(3/47r) 


Ss 


7.8125310"^ 


5.45402 10""^ 


5.34656 10"^ 


5.4016710"'^ 


7.9410"^ 


0.02366 


0.005040 


Sd 


7.8125310"^ 


5.0142610""^ 


4.8129910^^ 


4.8552310"*^ 


7.2210"^ 


0.02329 


0.005022 


Sh 


7.8125310"^ 


4.94734 10-"^ 


4.63675 10"^ 


4.64246 10"*^ 


7.0410"^ 


0.02321 


0.005019 



Table 1: The first four non zero Taylor coefficients of sd{x), ss{x) (a complex scalar) and sh{x) for 

X around vr, s{x) ~ Cj'^'' ~ '''")"' ) ^^'^ the coefficient of the term 1/x of these functions for x — > 0, 

s{x) ~ ci^l/x. The value of the functions for x = 7t/2 and x = 3/4 vr are also shown. sh{x) is normalized 
such that the quadratic coefficient coincide with the ones of soix) and ss{x). 
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Figure 2: Comparison between the coefficients of the logarithmic term for different theories. The picture 
shows, from top to bottom, the functions ss, sd and sh normalized to sd- The differences increase toward 
the origin where these ratios get ss/sd = 1-099 and sh/sd = -975 (not shown). 



where C [t) is the entropic c- function corresponding to a Dirac field and complex scalar in two dimensions 
respectively The integral in (j64p is done numerically using the analytic expressions for C{t) given in 

In figure 2 we have plotted ss{x) / sd{x) and sh{x) / sd{x) for a range x where our integration of these 
functions is precise enough. Clearly, they are linearly independent. This points to the effect of interactions 
in the CFT behind sh{x). 

As a cross check of these results one can verify in each case that s{x) satisfies the constraints coming 
from the strong subadditive inequality for the entropy. These are given by 

s{x) > 0, (65) 
s'{x) < 0, (66) 
s'{x) +sm{x)s"{x) > 0, (67) 

where x G (0, tt]. They are obtained by considering the strong subadditive inequality for different spatial 
(but relatively boosted to each other [6j ) plane angular sectors with a common vertex in Minkowski space. 
This problem can be mapped to a two dimensional de-Sitter space by taking the intersection of the plane 
angular sectors with an hyperboloid x'^x^ = — 1. Then, the inequalities (|65ti67p for the two dimensional 
de-Sitter case follow in a similar way as in the flat space case discussed in j5j (see also ^2j where a weaker 
form of ([67|) was derived). 

The coefficient of log(e) does not depend on the mass. Therefore we conclude that there are conformal 
field theories in three dimensions with non proportional entropy functions. This is in contrast to the 1 + 1 
case, where all the CFT would have entanglement entropy functions which are proportional to a unique 
function of the set V jl3J [23] (see however the recent works in [23] , and [5] ) . 

However, note that all the functions are remarkably similar to each other. There is a maximal relative 
difference between the scalar and Dirac case of 9% (the maximum relative distance seems to happen 
for X — > 0). The Hirata-Takayanagi formula approaches better the fermion than the boson case, with a 
maximum relative difference with respect to So (at the origin) which is remarkably small: only 2.5%. 
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Appendix A: Gap in the spectrum of the Dirac operator on the sphere 
with boundary conditions specified on a segment of a great circle 

In Section 2 we used that the operator Y D with boundary conditions on a segment of a great circle given 
by (I40p has no eigenvalues u in the interval (1/2, 3/2). In this Appendix we provide the proof of this fact. 
Since the spectrum is symmetric around i/ = 1 we have to show there are no eigenvalues in (1/2, 1]. First 
we note that due to ([^T|) the eigenvectors of corresponding to eigenvalue v satisfy 

/^^y = -v{v-l)^y. (68) 



For a = the functions are regular and (I68p implies due to the positivity of —A for regular functions 
that V ^ (0, 1). The case = 1 is also not possible since in this case the spinor only contains eigenvectors 
of A with eigenvalues and thus it is a constant spinor, giving v = Q instead. 

Now, as we increase a from to 1/2 the eigenvalues start moving. In order that at some point an 
eigenvalue enters to (1/2, 1] it must cross the point v = 1/2. At = 1/2 we have 

A$i/2 = ^«'i/2, (69) 

and the spinor diverges at the extreme points of the cut as d"^/^ (see eq. (jl4p ). This means, according 
to [H] (see also Section 3 in this paper), that it is possible to write 

where the rriij are unknown constant coefficients and the Si{6,(p) and Si{R{0,ip)) are the functions 
introduced in [11] and described also in Section 3. These must satisfy the eigenvalue equation (|69p and 
thus correspond to mass M = 1/2. If = 1/2 is eigenvalue of Y D the equation (|7Up has to hold for some 
coefficients, since adequately choosing them we can equate the divergent contributions on both sides of 
the equation. Then the difference between the left and right hand sides satisfies the Helmholtz equation 
(A - M2) = 0, and is singularity free, what means it is identically zero. 

Inserting (|70p in the eigenvalue equation for Y'Di and using the relevant differential equations for 
Si{9,{p) and Si{R{9,(p)) (see [H]) we arrive at 

+ (2a- 1)6 sin(x/2) = 0. (71) 

Here u, Pi and b are functions of a and x which satisfy the differential equations given in appendix B (for 
M = 1/2). With these one readily shows that there are no solutions for a E (0, 1/2). 

Appendix B: The Green function of a scalar field in a sphere with 
boundary conditions specified on a segment of a great circle 



In order to make this paper self contained we include the results of [TT] on the Green function of a complex 
scalar field on a sphere with the boundary conditions (P7|) . The trace is 

tr Gs{x, M, a) = 8^(1 - a)a [ Ha{y)dy . (72) 

J X 
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The function Ha{x) is the solution of the following set of ordinary non linear differential equations (we 
omit the subscript a and the dependence on x and the mass M of the variables for notational convenience) 

(73) 

(74) 
(75) 
(76) 
(77) 

(78) 

where Bi, B2, -B12, Pi, P2 are functions of x given in terms of H, Xi, X2, c, b, and u by the following set 
of algebraic equations 



H' = 


M 

-— {hB2 + cBi + 2uBi2) 






X[ = 


-M (bBu + uBi) , 






X'2 = 


-M {CB12 + UB2) , 






c = 


—2Mf32U csc(2;) sm(2;/2) — 


c (1 — a) 


csc(xj (1 + cos(xjj 


h' = 


-2M Piu csc{x) sin(2;/2) - 


ba csc(x 


) (1 + cos(x)) , 


u' = 


M 

-— sec(x/2) (6/32 + c/3i) + 


— u tan 1 
2 ^ 


x/2) , 



cos(x/2) 
87ra(l — a) 

sin(3;/2) 
87ra(l — a) 





sin(x/2) H - M {(3iX2+ P2X1) + 2M cos (x/2) u Bu , 

-cos(x/2) H -M tan(x/2) {P1X2 + P2X1) + M sm{x/2){bB2 + cBi) , 



(79) 
(80) 



-Msin(x/2)(cXi - 6X2) + Mtan(3;/2)(/32Si - P1B2) + (1 - 2a) cos(x/2)5i2 , (81) 
-4a(a - 1) - M2(4 - 8/3i/32 + bc + 3u^) 

-4 cos(x) (a(a - 1) +M2(ti2 + 1)) + cos(2x)(6c - n^) , (82) 
(2a- l)ticos(x/2) + Mtan(x/2)(/?ic-6/32). (83) 



The boundary conditions at x ^ vr are 



H{tt 

^2(vr 
u{tt 

bin 
cfvr 



0, 



r(-a) (cosh (^) Im U + a + 



2 



fsinh(^; 



22 '^/ii (cos (2a7r) + cosh(7r/ii)) r(l + a) 

Xi{tt) I a^(l-a) ) 

0, 



2 ' 



il-2a 



a(l — a) 



Mr2(l + a) 

^(■71") I a^(l-a) , 



(84) 
(85) 

(86) 
(87) 

(88) 
(89) 



where /ii = V 4 — 1 and is the digamma function. The meaning of the extra variables Bi, B2, B12, 
Xi, X2, n, 6, c, Pi and /32 is the same as in |15]. The trace in (|72]) is regularized such that it vanishes 
when X = vr, where there is no vertex point and no logarithmic term is present in the entropies. 
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